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1.  INTRODUCTION 

The  topic  of  nonlinear  acoustics  has  been  included  at  previous 
NATO  study  institutes  going  back  to  1966.  The  first  treatment,  by 
Berktay  (1),  speculated  on  the  possibilities  nonlinear  effects 
might  offer.  In  1968,  Mellen  (2)  illustrated  some  of  these  possi¬ 
bilities  through  laboratory  tank  experiments.  Berktay  returned  in 
1972  (3)  to  present  some  engineering  models  for  the  design  of  para¬ 
metric  sources.  In  1976,  Bj^rn^  (4)  presented  a  survey  of  theo¬ 
retical  and  experimental  results  on  parametric  arrays  developed  at 
several  laboratories. 

Today,  research  and  development  in  nonlinear  acoustics  has 
gravitated  toward  applications.  The  present  paper  therefore 
addresses  applications  with  a  view  towards  outlining  the  unique 
features  of  nonlinear  arrays,  especially  with  regard  to  signal 
processing.  Both  nonlinear  sources  and  receivers  will  be  discussed. 

As  research  tools,  nonlinear  sources  offer  a  means  of 
quantifying  both  the  limitations  of  the  medium  as  well  as  the 
payoffs  of  new  approaches  and  techniques.  In  the  following  pages, 
we  will  illustrate  this  claim  through  presentation  of  results 
obtained  with  several  nonlinear  research  tools  configured  for  model 
tests  as  well  as  full-scale  experimentation. 

Parametric  receivers  also  have  potential  as  research  tools 
but  they  are  of  interest  primarily  for  their  application  as  both 
fixed  and  mobile  sensors,  pending  the  successful  development  of 
practical  signal  processing  techniques. 
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2.  NONLINEAR  SOURCES 


As  is  well  known,  nonlinear  sources  create  highly  directive 
sound  through  the  nonlinear  interaction  of  intense  radiations  that 
are  themselves  directive.  Parametric  interaction  of  two  high  fre¬ 
quency  waves  yields  difference  frequency  waves,  while  the  dis¬ 
tortion  and/or  self-interaction  of  a  single  frequency  wave  leads 
to  the  creation  of  harmonic  radiations.  Difference  frequency 
radiations  have  the  advantage  of  low  frequency,  high  directivity, 
and  wide  bandwidth,  all  realized  from  a  relatively  small  primary 
source.  There  is  a  disadvantage  of  low  conversion  efficiency. 

Much  of  the  conversion  loss  from  primary  to  difference  frequency 
sound  can  be  made  up  by  signal  processing  gains,  as  will  be 
subsequently  discussed.  Harmonic  radiations  have  the  advantage 
of  high  directivity,  wide  bandwidth,  and  good  conversion  efficiency, 
with  the  disadvantage  of  high  frequency  that  is  prone  to  absorp¬ 
tion  by  the  medium. 

Those  not  familiar  with  nonlinear  acoustics  may  find  one  of 
the  tutorial  or  review  articles  useful  (4,5). 

2.1.  Directivity  of  Parametric  Sources 

The  beamwidth  capability  of  the  parametric  array  has  as  a 
theoretical  limit  the  Westervelt  result  (6), 


where  0^p  is  the  half-power  beamwidth,  ci  is  the  mean  absorption 
coefficient  of  the  primaries,  and  ks  is  the  difference  frequency 
wave  number.  Here  we  see  that  the  lower  the  absorption,  the 
longer  the  array  and  the  greater  the  directivity  of  the 
parametric  radiation. 

Eq.  (1)  is  shown  graphically  in  Fig.  1  for  a  downshift  ratio 
of  ten  (primary  frequency  f  *  difference  frequency  fs).  Actu¬ 
ally,  the  directivity  function  of  a  parametric  array  is  the 
convolution  of  the  product  of  the  primary  directivity  functions, 
Dp(0),  with  the  Westervelt  end-fire  array  directivity  (5),  Dw(0), 
weighted  by  sin<J),  i.e., 

Ds(0)  =  f  D  (<J>)  nw(d-<f>)  sin<J>d<J>  ,  (2) 

where  0  is  the  observation  angle  measured  from  the  acoustic  axis, 
4>  is  the  angle  variable  (to  an  annulus  of  integration  about  the 
acoustic  axis),  and 
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Figure  1.  Farfield  beamwidth  capability  for  a  para¬ 
metric  array  having  a  downshift  ratio  of  10,  as 
predicted  from  the  Westervelt  formula,  using  seawater 
absorption  data  as  indicated. 

Dw(0-(f>)  =  i  5  +  ks  sin2(^-=-^)  .  (3) 

Since  the  convolution  effect  accounts  for  broadening  of  the 
parametric  directivity  due  to  the  finite  size  of  the  primary 
beams,  the  Westervelt  directivity,  Eq.  (1),  and  the  corresponding 
plot  of  Fig.  (1)  can  only  be  approached  as  a  limit.  It  can 
nonetheless  be  seen  that  warm  water  permits  smaller  beamwidths 
than  cold  water  (Shulkin  and  Marsh  absorption  data)  while  con¬ 
sideration  of  the  boron-borate  relaxation  process  (Thorp  data)  (7) 
leads  to  practically  constant  parametric  directivity  below  about 
5  kHz.  If  maximum  parametric  directivity  is  desired,  primary 
frequencies  near  5  kHz  appear  optimum. 

2.2.  Applications  of  Directive  Parametric  Devices 

Possession  of  a  narrow  system  beamwidth  with  good  minor  lobe 
suppression  enables  one  to  reject  environmental  backscatter  and 
improve  echo-to-reverberation  ratios,  E/R.  Usually,  improved  E/R 
increases  detection  probability,  although  just  being  able  to 
detect  the  target  often  provides  classification  clues,  as  is  the 
case  in  Doppler  measurements.  Extremely  narrow  beamwidths  are 
useful  in  acoustic  imaging,  which  enables  the  classification  of 
targets  as  to  size  and  shape. 

We  argue  here  that  these  capabilities  are  equivalent  to 
signal  processing  capabilities  since  the  same  goals  are  achieved 
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in  both  cases.  We  therefore  have  signal  processing  upon 
transmission  with  a  directive  nonlinear  source. 

2.2.1.  Model  submarine  experiments.  An  example  of  the 
utility  of  the  high  directivity  of  the  parametric  array  in  sonar 
research  is  offered  by  the  data  of  Fig.  2.  Here  we  see  some 
normalized  scattering  patterns  from  a  model  submarine  some  2  m 
in  length  that  were  acquired  with  the  target  placed  in  rotation 
at  mid  depth  in  a  shallow  water  lagoon.  The  water  was  1.4  m  deep, 
overlaying  a  sand  bottom.  The  range  was  250  m,  and  the  signal 
was  a  3.5  msec  cw  ping  at  15  kHz.  Two  patterns  are  shown.  One 
is  for  a  linear  radiation  from  an  array  23  cm  in  diameter  with  a 
half-power  beamwidth  of  26°.  The  other  was  acquired  with  a 
parametric  array  developed  from  primaries  centered  at  225  kHz  and 
emitted  through  the  same  sized  aperture  to  develop  a  parametric 
beam  3°  wide.  As  expected,  the  broad  beam  linear  system  insoni- 
fies  much  more  of  the  backscattering  surfaces  and  volume  scat- 
terers,  to  the  extent  that  the  target  is  really  only  detected 
near  beam  aspect.  The  narrow  beam  parametric  system  insonifies 
fewer  scatterers  and  the  mean  reverberation  is  suppressed  by 
11.9  dB  over  that  of  the  linear  system.  Simple  theory,  based 
solely  on  insonified  area,  would  predict  a  suppression  of  8.2  dB. 
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Figure  3.  Signal  envelope  for  model  submarine  echo  data. 


As  can  be  seen,  the  3°  parametric  beam  has  sufficient  E/R  to 
detect  the  target  at  practically  all  aspects. 

It  is  instructive  to  use  the  experimental  data  to  determine 
the  detection  probability  of  the  model  target  as  a  function  of 
detection  threshold.  The  latter  is  defined  for  a  50%  false  alarm 
rate  as 


D 


where 
T  = 


v*(t) 


dt 


(4) 


the  target  plus  reverberation  energy  in 
each  echo. 


V  (t)  dt>  .. 
r  ensemble 


the  reverberation  only  signal, 
and 


a 


2 

B 


dt>, 


ensemble 


the  variance  of  reverbera¬ 
tion  energy. 


where  Vr(t)  is  the  received  signal  and  the  times  t^  are  shown  in 
Fig.  3.  One  computes  B  and  as  constants  from  the  data  ensem¬ 
ble  and  then  computes  the  probability  of  observing  the  target 
as  the  ratio  of  number  of  times  each  measured  D  exceeds  a  given 
Dx  to  the  number  of  estimates  in  the  ensemble,  i.e., 

nD  >D 

Probability  [D^D  ]  =  — 5  .  (5) 


The  narrowbeam  parametric  system  achieved  a  95%  detection  prob¬ 
ability  at  a  0  dB  threshold,  compared  to  50%  for  the  broadbeam 
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Figure  4.  Detection  probability  for  model  submarine 
echo  data. 


case;  the  results  are  shown  in  Fig.  4.  With  increase  in  threshold, 
the  detection  probability  of  the  narrowbeam  system  falls  off 
slower  than  for  its  broadbeam  counterpart. 

2.2.2.  Doppler  measurements.  Another  example  of  the  value 
of  parametric  directivity  is  afforded  by  application  of  the  same 
apparatus  to  Doppler  detection.  Since  most  everything  of  interest 
in  the  ocean  is  in  movement,  we  are  obliged  to  understand  the 
dynamics  of  the  process.  With  respect  to  the  medium,  this 
involves  relating  the  directional  spectra  of  scatterers  in  motion 
to  the  directional  properties  and  motion  of  the  sonar  in  order  to 
arrive  at  the  Doppler  spectra  it  will  observe.  Although  this  pro¬ 
cess  is  intuitively  understood,  it  is  by  no  means  delineated,  and 
important  work  along  these  lines  is  in  progress  (8). 

Conceptually,  backscattered  Doppler  received  at  some  angle  6 
on  a  sonar  beam  can  be  modeled  by 

Af  =  2 ( f  /c  ) (v  -v  )  •  D(0)e  ,  (6) 

o  o  s  m  o 

-V 

where  fQ  is  the  carrier  frequency,  cQ  is  the  sound  speed,  vs  is 
the  velocity  vector  of  the  ship  carrying  the  sonar,  vm  is  the 
velocity  vector  of  the  segment  of  the  medium  in  view  (which  may 
contain  a  target),  D(9)  is  the  directivity  function  of  the  sonar, 
and  <£  is  the  unit  vector  of  the  angular  position  on  the  beam 
relative  to  that  of  the  ship's  path.  This  simplified  result  would 
give  the  spectrum  of  backscattered  Doppler,  if  all  of  the 
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scatterers  in  the  medium  had  the  same  target  strength.  Actually, 
vm  is  itself  a  spectrum  that  can  be  quite  complicated  for  the  case 
of  a  turbulent  medium.  Since  each  inhomogeneity  there  has  its  own 
cross  section  as  well  as  velocity,  rigorous  solutions  of  the  prob¬ 
lem  are  usually  expressed  in  terms  of  a  scattering  integral. 

Eq.  (6)  nonetheless  illustrates  an  important  aspect  of 
Doppler  measurements  that  involves  the  influence  of  the  direc¬ 
tivity  function  D(0).  This  function  weights  Doppler  frequency 
measurements  in  reverberant  environments  such  that  the  broader 
the  beam  the  broader  the  Doppler  spectrum  received.  In  this 
case,  precise  Doppler  data  can  only  be  acquired  with  systems  that 
have  good  directivity. 

Some  experimental  results  on  this  topic  are  depicted  in 
Fig.  5.  A  spherical  target  was  towed  at  mid  depth  toward  a  fixed 
soundhead  operating  at  20  kHz  in  the  aforementioned  lagoon  experi¬ 
ment.  The  target  was  moving  at  approximately  2  kt.  The  surface 
of  the  water  was  driven  by  the  wind  at  a  comparable  velocity 
along  the  same  direction  of  motion.  A  50  msec  cw  pulse  at  15  kHz 
was  used  to  insonify  the  target  for  Doppler  measurements.  In 
processing,  the  received  signal  was  beat  against  the  transmit 
signal,  and  the  Doppler  shifted  frequencies  were  obtained  at 
baseband  by  low  pass  filtering.  These  are  displayed  in  the  figure 
in  quasi-3D  plots  that  show  the  range  dependence  of  Doppler 
signal  amplitude  (on  the  vertical  axis)  for  successive  pings 
(time  history)  as  the  target  closes  range. 

It  can  be  seen  that  the  target  appears  clearly  in  the 
narrowbeam  data.  Its  Doppler  frequency  is  also  measurable  from 
the  baseband  oscillations  along  the  range,  or  equivalently,  the 
instantaneous  time  axis.  The  broadbeam  system  views  more  of  the 
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Figure  5.  Doppler  data  from  a  moving  target  masked  by 
moving  surface  reverberation. 
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moving  surface  and  therefore  experiences  more  Doppler  noise  which 
masks  the  target  and  makes  detection  and  measurement  more  difficult. 

2.3.  Directivity  of  Harmonic  Radiations 

Another  important  directivity  aspect  of  nonlinear  acoustics 
lies  in  the  area  of  harmonic  radiations  created  in  the  nonlinear 
distortion  and  shock  formation  process.  As  is  well  known,  the 
velocity  varies  within  a  high  intensity  waveform  in  proportion  to 
the  local  particle  velocity.  This  causes  a  steepening  of  the 
waveform  in  the  time  domain,  which  is  equivalent  to  the  creation 
of  harmonic  components  in  the  frequency  domain. 

The  harmonics  are  attractive  because  of  their  great 
directivity.  It  has  been  shown  (9)  that  the  directivity  function 
of  the  nth  harmonic  is  equal  to  that  of  the  fundamental  raised  to 
the  nth  power,  or 


D  (0)  =  D"(0)  .  (7) 

n  1 

Thus,  as  one  progresses  to  the  higher  harmonics,  the  width  of  the 
beam  decreases  and  the  minor  lobe  suppression  increases.  These 
are  desirable  features  in  imaging  applications. 

The  amplitudes  of  the  harmonic  radiations  can  be  quite  high. 
At  initial  shock  formation,  the  2nd,  3rd,  and  4th  harmonics  are 
only  8,  12,  and  15  dB,  respectively,  below  the  fundamental.  In 
hard  shock,  their  amplitudes  go  as  the  reciprocal  harmonic  number, 
so  that  the  same  respective  harmonic  amplitudes  will  be  6,  10, 
and  12  dB  below  the  fundamental  (9). 

2.4.  Applications  of  Directive  Harmonics 

We  have  conducted  several  experiments  with  harmonic  sonar 
using  a  rectangular  transducer  measuring  1.0  x  0.1  m.  It  was 
energized  at  a  fundamental  frequency  of  100  kHz  in  200  psec  cw 
pulses  with  1.3  kW  of  pulse  power.  With  a  transducer  efficiency 
of  50%,  this  yields  a  source  intensity  of  0.7  W/cm^.  Since  the 
beam  is  only  0.8°  x  8°  in  half-power  width,  it  has  a  high 
directivity  and  source  level,  which  means  that  harmonic  radiations 
are  created  as  the  high  intensity  pulse  propagates  through  the 
water.  Echoes  from  targets  are  received  by  a  wideband  hydrophone, 
amplified,  and  passed  through  any  of  a  series  of  narrow  bandpass 
filters  tuned  to  the  harmonic  frequencies.  The  harmonic  signals 
are  processed  and  displayed  in  a  range/bearing  (B  scan)  format. 

2.4.1.  Detection  sonar.  A  typical  illustration  of  some  of 
the  advantages  of  harmonic  echo  scanning  is  shown  in  Fig.  6.  A 
point  target  at  a  range  of  100  m  is  depicted  with  both  the 
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Figure  6.  Comparison  of  point  target  detection 
performance . 


fundamental  and  its  second  harmonic  component.  The  fundamental 
radiation  has  a  half-power  resolution  of  0.8°  or  1.4  m  in  bear¬ 
ing  and  0.3  m  in  range,  with  a  minor  lobe  suppression  of  13  dB. 

A  scan  taken  at  the  fundamental  frequency  shows  the  target  in  a 
condition  of  display  saturation,  which  is  a  quite  common  occur¬ 
rence  in  sonar  displays.  The  second  harmonic  data  demonstrate 
the  alleviation  of  these  undesirable  effects,  even  though  the 
video  gain  was  normalized  to  the  same  peak  input  voltage.  With 
the  second  harmonic,  the  half-power  beamwidth  was  reduced  tc  0.3°, 
yielding  a  resolution  cell  0.9  m  in  azimuth  by  0.3  m  in  range. 

Also,  the  minor  lobe  suppression  was  increased  to  around  20  dB 
(it  should  have  been  26  dB).  The  target  display  is  much  more 
compatible  with  the  sonar's  resolution  capability  and  the  adjust¬ 
ment  of  the  display  is  much  less  critical. 

2.4.2  Imaging  Sonar.  Further  demonstration  of  the  dynamics 
of  harmonic  echo  scanning  is  afforded  by  the  data  of  Fig.  7  which 
depicts  the  viewing  of  larger  targets.  Here  the  beam  patterns  are 
shown  at  left  for  the  1st  through  5th  harmonics  of  the  experiment 
at  hand.  Range/bearing  images  of  a  cylinder  are  shown  at  right, 
with  each  display  corresponding  to  a  particular  harmonic  echo  scan, 
each  made  at  a  range  of  100  m.  With  the  fundamental  radiation,  the 
beam  is  almost  as  large  as  the  target.  As  a  result,  one  really 
sees  only  some  large  "blobs",  probably  due  to  echoes  from  each 
end  of  the  target.  Reception  of  these  echoes  with  the  minor  lobes 
is  also  evident.  The  sidelobe  echoes  are  not  a  problem  in  the 
second  harmonic  data.  With  the  third  harmonic,  the  orientation 
of  the  target  becomes  clearly  evident.  Data  for  the  fourth 
harmonic  shows  a  fair  representation  of  the  size  of  the  target. 

This  trend  continues  to  the  fifth  harmonic,  although  the  display 
has  become  somewhat  sparse,  perhaps  for  nonacoustic  reasons. 

The  data  displayed  in  the  center  column  of  Fig.  7  are  for 
the  viewing  of  an  object  much  larger  than  the  original  beam. 
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200  kHz 
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BEAM  PATTERNS  FLOATING  BARGE  SUBMERGED  CYLINDER 

Figure  7.  Imaging  with  the  harmonic  radiations. 


These  are  range/bearing  scans  of  a  barge  having  a  diving  plat¬ 
form  with  pilings  attached.  As  can  be  seen,  the  definition  in 
the  field  of  view  becomes  progressively  better  with  increase  in 
harmonic  number,  angular  resolution,  and  minor  lobe  suppression. 

At  about  the  fourth  harmonic,  the  four  pilings  (15  cm  in  diameter) 
becomes  clearly  evident. 


Although  it  can  be  argued  that  the  harmonic  radiations  are 
subject  to  more  absorption  than  the  fundamental,  it  is  nonethele.^ 


true  that  many  applications  are  not  precluded  by  absorption 
limited  effects.  It  appears  reasonable,  for  example,  to  acquire 
crude  images  with  the  low  frequency  fundamental  for  initial  study 
and  then  develop  sharper  images  with  the  harmonics,  as  the  range- 
dependent  absorption  factor  permits. 

2.5.  Bandwidth  of  Parametric  Sources 

The  frequency  response  of  a  parametric  array  may  be  obtained 
from  the  principle  of  superposition  applied  to  any  of  the  steady 
state  solutions  for  difference  frequency  sound.  This  has  been 
formulated  (10)  for  the  Westervelt  solution  (6)  on  the  acoustic 
axis  as 


CO 

p  (w  )  11  to  /  P,  (to)  p„(w+w  )dco 

S  S  S  #  J.  Z  S 

•'  —  CO 


where  ws  is  the  difference  frequency  and  Pj(')  and  P2O)  ate  the 
two  primary  (interacting)  pressures.  This  result  is  equivalent 
to  a  crosscorrelation  integral  in  product  with  the  difference 
frequency  squared.  It  is  a  very  straightforward  and  powerful  re¬ 
sult  that  has  been  used  with  the  theorems  of  convolution,  Fourier 
transform,  etc.,  to  obtain  the  transient  response  of  the  parametric 
array  (10). 

This  result  can  also  be  used  to  determine  the  difference 
frequency  response  of  a  given  system  by  considering  the  frequency 
response  of  the  primary  radiations  in  their  own  operating  band. 
Primary  projectors  are  almost  invariably  tuned  devices  that  are 
impedance  matched  to  their  power  amplifiers  at  their  resonance 
frequencies.  Each  projector  has  a  limited  Fourier  spectrum  over 
which  it  can  transmit  sound  at  high  intensity.  Consideration  of 
the  steady  state  frequency  response  enables  Eq.  (8)  to  be  reduced  to 

dp  Ju>  ) 

— ^ pi(Ui)  ,  (9) 


where  Ws=u)l-W2 ’  anc*  ul>  w2  denote  the  primary  frequencies  utilized. 
To  treat  a  practical  example,  let  us  suppose  that  we  have  a  typical 
projector  in  the  primary  frequency  band  whose  half-power  bandwidth, 
A«hp>  is  about  15%  of  the  resonance  frequency,  u)Q,  and  whose 
response  is  a  Gaussian  function  of  frequency,  i.e., 


2  2 

=  p  exp  -a  (to-w  ) 
o  o 


where  a  =  7.85  sec/rad,  a  constant.  A  projector  of  this  type 
would  have  a  Q=w  /Aw  _  of  6.7.  This  function  is  plotted  at  left 
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Figure  8.  Parametric  bandwidth  and  frequency  response 
characteristics  for  a  typical  primary  source. 


in  Fig.  8.  We  further  assume  that  the  two  primary  radiations 

are  driven  symmetrically  about  the  resonance  frequency  u>o,  as  shown. 

Application  of  Eq.  (9)  to  the  difference  frequency  response 
yields  the  plot  at  right  in  Fig.  8.  It  can  be  seen  that  the  dif¬ 
ference  frequency  pressure  maximizes  at  £>s%18%  of  the  primary 
frequency.  There  is  essentially  a  bandwidth  translation  from  pri¬ 
mary  to  difference  frequency  sound,  as  shown.  In  terms  of  ws, 
the  difference  frequency  radiation  would  be  capable  of  about  an 
85%  bandwidth,  or  a  Q  of  0.2.  These  features,  as  well  as  the 
asymmetric  shape  of  the  difference  frequency  spectrum,  are  a  re¬ 
sult  of  the  dependence  of  Eq.  (8)  on  the  p^(*)  P2(*)  product,  up¬ 
ward  weighted  by  multiplication  with  As  a  result,  the 

difference  frequency  is  capable  of  a  -3  dB  bandwidth  greater  than 
one  octave,  and  is  useful  over  even  larger  bandwidths. 

2.6.  Processing  of  Wideband  Parametric  Signals 

Given  that  a  parametric  sonar  will  usually  have  greater 
directivity  and  bandwidth  than  a  comparable  linear  device,  it  is 
worthwhile  to  consider  the  difference  this  may  make  in  signal 
processing.  These  attributes  may  be  utilized  in  various  ways, 
from  measuring  the  frequency  response  of  a  resonant  target  in  a 
reverberant  environment,  to  extracting  a  weak  signal  from  noise 
or  reverberation  with  wideband  processing,  etc.  Since  it  is  not 
our  purpose  here  to  invent  new  active  sonar  processing  techniques. 
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Figure  9.  Spectra  of  a  wideband  FM  chirp  at  several 
ranges  in  a  shallow  water  waveguide. 


we  will  restrict  the  discussion  to  existing,  well-known 
procedures. 


2.6.1.  Spectral  measurements.  The  bandwidth  capability  of 
a  parametric  array  may  be  used  to  obtain  the  frequency  response  of 
a  wide  variety  of  resonant  structures  as  well  as  the  spectral 
response  of  environmental  features  (reverberation,  model  propaga¬ 
tion,  etc.).  The  FM  chirp  signal  is  especially  suited  to  rapid, 
and  practically  instantaneous,  acquisition  of  this  type  of  data. 

An  example  is  afforded  by  the  results  of  Fig.  9,  which  shows 
some  normalized  propagation  spectra  associated  with  the  trans¬ 
mission  of  a  wideband  parametric  pulse  in  the  previously  described 
lagoon  experiment.  Here  the  source  and  receiver  were  at  mid  depth 
in  the  1.4  m  waveguide  and  the  signal  was  swept  over  a  little  more 
than  an  octave  in  an  FM  chirp  3.5  msec  long.  The  wind  driven  sur¬ 
face  waves  during  the  acquisition  of  data  had  a  standard  deviation 
2a  =  6.2  cm.  Spectra  measured  with  an  analyzer  bandwidth  of  1  Hz 
are  shown  at  three  ranges.  At  27  m,  the  spectrum  has  a  half¬ 
power  width  of  21  kHz,  but  it  falls  to  18.5  kHz  at  320  m  and  16  kHz 
at  640  m,  which  translates  into  a  loss  in  bandwidth  of  8.2  Hz  per 
meter  of  propagation  range.  Jensen  and  Kuperman  (11)  have  sug¬ 
gested  that  the  rounding  off  of  the  spectrum  at  the  lower  fre¬ 
quencies  is  due  to  bottom  loss,  while  losses  at  the  upper  edge  of 
the  spectrum  may  be  due  to  volume  and  surface  scatterers.  The 
capacity  of  the  medium  to  sustain  wide  bandwidth  signals  is  an 
important  issue  in  underwater  acoustics.  It  is  interesting  and 
useful  to  study  the  governing  mechanisms  with  different  research 
tools . 
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Figure  10.  Implications  for  processing  gain  as  a 
function  of  attainable  time-bandwidth  properties. 


2.6.2.  Large  time-bandwidth  processing.  For  the  simple 
processors  in  common  use,  one  can  talk  in  terms  of  coherent  or 
incoherent  techniques.  The  latter,  in  the  form  of  linear  or  square 
law  detector/averager  systems,  may  be  preferred  for  distributed 
targets  or  where  the  medium  multipath  destroys  the  coherence  of 
the  signal.  The  former,  in  the  form  of  matched  filter  systems 
(replica  correlator,  delay  line,  or  other  implementation)  may  be 
preferred  when  the  degradation  of  signal  coherence  is  not  pro¬ 
hibitive.  Both  approaches  have  in  common  the  fact  that  the  pro¬ 
cessing  gain  is  proportional  to  the  time-bandwidth  product,  TW, 
of  the  signal. 

To  illustrate  the  capability  of  parametric  arrays  in  this 
regard,  we  assume  a  coherent  processor  whose  processing  gain 
against  white  Gaussian  noise  is  (12) 

P.G.  =  10  log  TW  .  (11) 


Some  results  for  an  assumed  center  frequency  of  3  kHz  are  shown 
in  Fig.  10  as  a  function  of  pinglength  for  15%,  1/3  octave,  and 
1  octave  bandwidths.  There  is  a  3.5  dB  improvement  in  going  from 


15%  bandwidth  to  1/3  octave,  and  an  8.2  dB  improvement  in  going 
to  a  1  octave  bandwidth.  Whether  or  not  the  medium  or  the  target 
can  accommodate  large  bandwidths  is  an  issue  worthy  of  resolving. 
If  such  large  bandwidths  are  attainable,  they  will  probably  occur 
at  low  frequencies,  where  propagation  losses  are  minimized  and 
where  low  ka  target  strengths  can  be  realized.  Here,  the  targets 
would  behave  more  like  point  targets  with  fairly  omnidirectional 
scattering  patterns.  Although  it  is  difficult  to  realize  low  ka 
systems,  those  with  ka=l-3  appear  attractive.  The  parametric 
array  approach  to  low  ka  systems  utilizes  nonlinear  acoustics' 
strongest  attribute — low  frequencies  with  small  transducers. 


If  the  noise  is  reverberation,  there  is  an  added  signal-to- 
noise  gain  that  is  equivalent  to  a  processing  gain.  This  can 
easily  be  included  in  Eq.  (11)  if  we  recognize  that  at  low  audio 
frequencies,  the  beamwidth  of  the  parametric  array  will  be  domi¬ 
nated  by  the  product  directivity  of  the  primaries,  Dp,  in  the 
convolution  of  Eq.  (2).  The  difference  frequency  beamwidth  will 
then  be  slightly  larger  than  the  product  pattern  beamwidth.  For 
surface  reverberation,  the  processing  gain  improvement  in  going 
from  a  conventional  linear  system  to  a  parametric  system  (of  the 
same  size  operating  at  the  same  frequency)  is  proportional  to  the 
ratio  of  their  respective  insonified  areas.  This  is 
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where  the  downshift  ratio  fp/fs  is  again  the  ratio  of  primary  to 
difference  frequencies.  The  dashed  lines  in  Fig.  10  represent 
the  relative  merit  of  parametric  directivity  for  this  case. 

Directivity,  of  course,  is  as  useful  as  bandwidth  for 
processing  gain  in  a  reverberant  situation  and  it  may  be  easier 
to  achieve  and  maintain.  More  importantly,  it  may  often  be  the 
essential  factor  that  precludes  all  other  considerations.  As 
Lavai  points  out  (13),  the  detection  of  a  signal  in  reverberation 
is  either  a  "go"  or  "no  go"  situation.  One  simply  has  to  have 
enough  suppression  to  recognize  the  target  in  the  first  place, 
irrespective  of  the  source  level  or  efficiency. 

3.  NONLINEAR  RECEIVERS 

The  same  general  principles  at  work  in  generation  of  low 
frequency  signals  by  nonlinear  sources  are  involved  in  reception 
of  low  frequency  acoustic  waves  by  a  nonlinear  receiver.  In  the 
parametric  source,  information  in  the  two  high  frequency  primary 
waves  is  translated  by  intermodulation  to  a  lower  frequency  region 
because  of  the  inherent  nonlinearity  of  the  pressure  density 
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relationship  of  water.  On  the  other  hand,  in  the  nonlinear 
receiver  the  information  in  a  low  frequency  signal  wave  is  trans¬ 
lated  to  the  sidebands  of  the  pump  wave  by  intermodulation  of  the 
low  frequency  signal  wave  and  the  pump  wave. 

In  the  parametric  source,  essentially  all  of  the  signal 
processing  is  performed  by  the  water  since  differential  absorption 
removes  the  high  frequency  primary  waves  and  leaves  only  the  low 
frequency  signal  that  propagates  to  long  ranges.  However,  in  the 
case  of  the  nonlinear  receiver,  the  desired  information  exists  as 
low  level  modulation  sidebands  of  the  pump  signal,  which  places 
severe  constraints  on  the  signal  processing  used  to  recover  the 
information.  Signal  processing  techniques  are  the  focus  of  the 
remainder  of  this  paper,  but  to  introduce  that  topic,  some  back¬ 
ground  material  about  the  nonlinear  receiver  will  be  given  first 
to  establish  a  common  base  of  terminology. 

3.1.  Background 

The  nonlinear  acoustic  receiving  array,  or  parametric 
receiver,  uses  the  inherent  nonlinearity  of  the  water  to  achieve 
directional  reception  of  low  frequency  acoustic  waves  with  only 
two  small  high  frequency  transducers  and  some  associated  elec¬ 
tronic  hardware.  This  concept  was  first  suggested  by  Westervelt, 
almost  as  an  afterthought,  in  his  classic  paper  on  parametric 
acoustic  arrays  (6).  In  the  next  few  years,  a  number  of  theo¬ 
retical  and  experimental  investigations  followed,  most  of  which 
emphasized  demonstrating  existence  of  the  phenomenon  and  de¬ 
veloping  and  validating  mathematical  models  to  describe  the  basic 
physics  of  the  process,  e.g.,  Truchard  (14). 

The  basic  elements  of  the  parametric  receiver  are  the  pump 
oscillator  and  power  amplifier,  two  high  frequency  transducers, 
and  the  receiving  electronics.  The  pump  oscillator  and  power 
amplifier  generate  a  stable,  continuous,  high  frequency  signal 
and  amplify  it  to  a  level  sufficient  to  produce  the  desired  pump 
acoustic  level  in  the  water.  The  pump  signal  is  projected  from 
one  of  the  transducers  (pump)  and  received  by  the  second  trans¬ 
ducer  (hydrophone)  located  a  distance  L  from  the  pump.  Ambient 
low  frequency  acoustic  waves  propagating  through  the  same  water 
volume  will  interact  nonlinearly  with  the  pump  wave  to  generate 
intermodulation  products.  The  function  of  the  receiver  electronic 
hardware  is  to  recover  the  information  contained  in  the  ambient 
signals  by  demodulating  the  interaction  products  which  appear  as 
modulation  sidebands  on  the  pump  carrier.  Thus  an  ambient  acous¬ 
tic  wave  of  frequency  fs  produces  an  electrical  signal  at  the 
output  of  the  receiver  electronics  which  is  also  of  frequency  fs. 

Zverev  and  Kalachev  (15)  expressed  the  pump  signal  and  the 
intermodulation  products  at  the  hydrophone  in  closed  form  as  a 
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phase  modulation  of  the  pump  wave  by  the  ambient  acoustic  waves. 
If  the  pump  wave  is  an  ideal  sinusoid  cos(wQt),  the  signal  at 
the  hydrophone  can  be  represented  as 

s(t)  =  cos  (a)  t+<J>  (t)+<jj)  ,  (13) 

o  a 

where  <j>a(t)  is  a  phase  modulation  due  to  the  nonlinearity  of  the 
water  and  iJj  is  a  phase  shift  related  to  the  propagation  delay 
from  the  pump  to  the  hydrophone.  The  acoustic  phase  modulation 
term  is  given  by 


<f>  (t) 
a 


[6  -  (l-cos0)]  Pg  sin[  (kgL/2)  (l-cos0)  ] 

3  (k  L/2) (1-cos0) 

poco  s 


(14) 


In  this  expression  0  is  the  plane  angle  measured  from  the  line 
joining  the  pump  and  hydrophone,  Pg  and  kg  are  the  amplitude  and 
wave  number,  respectively,  of  the  signal  to  be  detected,  L  is  the 
pump-hydrophone  separation,  w0  is  the  pump  frequency,  p0  and  c0 
are  the  equilibrium  density  and  sound  speed,  respectively,  and 
B(=1+B/2A)  is  the  coefficient  of  nonlinearity  of  the  medium 
(approximately  equal  to  3.5  in  seawater). 


The  directional  response  of  the  parametric  receiver  can  be 
obtained  from  Eq.  (14)  and  is  given  by 


D(0) 


6  -  (l-cos0)  sin[(ksL/2)(l-cos0)] 
3  (ksL/2) (1-COS0) 


(15) 


From  Eq.  (15)  it  is  clear  that  the  directional  response  of  the 
parametric  receiver  is  symmetric  about  the  maximum  response  axis 
which  is  in  the  direction  of  a  line  extending  from  the  hydrophone 
through  the  pump.  Furthermore,  the  directional  response  is 
independent  of  the  pump  frequency;  it  depends  only  upon  the  pump- 
hydrophone  separation  and  the  wavelength  of  the  signal  to  be 
detected.  Although  the  array  synthesized  in  the  interaction 
volume  is  actually  a  volumetric  array,  the  directivity  character¬ 
istics  of  the  synthesized  array  are  similar  to  those  of  a 
continuous,  end-fire  array  of  length  L.  Thus,  the  half-power 
beamwidth  is  given  in  radians  approximately  by 

0Hp  =  1.9iA8/L  ,  (16) 

where  As  is  the  acoustic  wavelength  of  the  signal  to  be  detected. 
It  is  the  end-fire  array  effect  that  provides  the  directivity 
of  the  parametric  receiver  and  hence  its  ability  to  discriminate 
against  the  low  frequency  ambient  noise  that  may  otherwise  mask 
the  signal  wave. 
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The  detection  of  low  frequency  signals  in  the  ocean  is 
closely  related  to  the  ability  of  the  acoustic  sensor  to  dis¬ 
criminate  against  low  frequency  ambient  noise  and  thereby  im¬ 
prove  the  signal-to-noise  ratio  (S/N)  compared  to  a  simple, 
omnidirectional  sensor.  One  measure  of  the  S/N  improvement  of 
an  acoustic  sensor  is  spatial  processing  gain  (SPG),  which  is  the 
ratio  of  the  noise  received  by  an  omnidirectional  sensor  to  the 
noise  received  by  the  directional  acoustic  sensor.  If  the  noise 
field  is  isotropic,  the  SPC  is  equal  to  the  directivity  index  (DI) 
of  the  sensor.  For  large  acoustic  apertures  the  DI  of  the  para¬ 
metric  receiver  asymptotically  approaches  the  DI  of  a  continuous 
end-fire  array  of  the  same  acoustic  aperture, 

DI  =  10  log(4L/As)  .  (17) 

Although  the  ambient  noise  field  is  rarely  isotropic,  the  DI  is  a 
convenient  and  useful  measure  for  first  order  comparisons  of 
different  acoustic  sensors. 

Since  the  ambient  noise  field  is  usually  anisotropic,  the 
response  of  the  sensor  to  noise  from  the  back  side  is  very  impor¬ 
tant.  For  L/Ag>>l,  the  ratio  of  the  maximum  response  of  the  para¬ 
metric  receiver  to  the  envelope  of  the  back  lobes  is  given  in 
decibels  by 

F/B  =  20  log(147rL/3As)  .  (18) 

Both  the  DI  and  the  F/B  of  the  parametric  receiver  are 
functions  of  the  acoustic  aperture  and  are  independent  of  the 
pump  frequency. 

Any  acoustic  sensor  will  have  some  self-noise  floor  which 
represents  the  noise  pressure  equivalent  to  the  internal,  or 
self-noise,  sources.  Clearly,  the  self-noise  of  a  sensor  must  be 
less  than  the  amplitude  of  the  signal  to  be  received  by  the  sensor. 
The  self-noise  contributors  in  the  parametric  receiver  must  be 
identified  and  controlled  if  a  S/N  improvement  commensurate  with 
the  directivity  characteristics  is  to  be  achieved.  The  principal 
noise  sources  in  the  parametric  receiver  have  been  identified  and 
selection  of  parameter  values  to  minimize  the  effects  of  these 
noise  sources  have  been  discussed  previously  (16-19). 

The  difficulty,  and  the  importance,  of  the  receiver  elec¬ 
tronics  design  can  be  appreciated  from  the  following.  Since  the 
water  is  only  weakly  nonlinear,  the  interaction  components  (mod¬ 
ulation  sidebands)  are  very  low  level  signals.  For  spherically 
spreading  pump  waves,  Berktay  and  Shooter  (20)  derive  an  ex¬ 
pression  for  the  axial  value  of  the  interaction  component  pressure 
amplitude  as 
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Figure  11.  Near-sideband  signal  spectrum  of  the 
parametric  receiver  showing  the  upper  and  lower  sideband 
signals. 


w  BP  P 

p(L)  =  — — —  exp  (-a  L)  ,  (19) 


where  P  is  the  pressure  amplitude  of  the  pump  wave  at  unit 
distance  from  the  pump,  aQ  is  the  absorption  coefficient  at  the 
pump  frequency  w0,  and  we  have  made  use  of  the  fact  that 
wo>>u>s  to  simplify  the  expression. 

Berktay  and  Muir  (21)  define  pump  excess  as  the  ratio  of  the 
amplitude  of  the  pump  signal  to  the  amplitude  of  the  interaction 
component,  measured  at  the  hydrophone.  This  is  illustrated  in 
Fig.  11  which  shows  a  hypothetical  spectrum  at  the  output  of  the 
hydrophone  when  the  parametric  receiver  is  receiving  an  acoustic 
wave  of  frequency  fs.  The  large  amplitude  component  at  frequency 
fp  represents  the  pump  signal.  The  interaction  components 
denoted  "USB"  and  "LSB"  represent  the  upper  and  lower  sidebands 
resulting  from  modulation  of  the  pump  wave  by  the  acoustic  wave 
of  frequency  fg.  For  cases  of  practical  interest,  the  amplitude 
of  the  interaction  components  will  be  comparable  to  the  low  fre¬ 
quency  ambient  noise  within  the  parametric  receiver  beam.  Under 
these  conditions,  the  pump  excess  of  a  parametric  receiver  for 
passive  sonar  application  to  receive  low  audio  frequencies  in  the 
ocean  will  be  in  the  140  to  170  dB  range.  Thus,  although  the 
nonlinear  acoustic  process  is  phase  modulation,  we  can  ignore  the 
higher  order  sideband  components,  approximate  the  process  as  a 
linear  modulation,  and  use  linear  modulation  techniques  to  recover 
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Figure  12.  Block  diagram  of  the  bandpass  receiver. 


the  sideband  signals.  The  modulation  index  is  then  equivalent 
to  the  reciprocal  of  the  pump  excess,  and  it  is  clear  that  a 
receiver  capable  of  recovering  sideband  signals  with  modulation 
indices  in  the  range  of  10-7  to  3*10-9  will  be  required. 


3.2.  Sideband  Recovery  Techniques 


Three  techniques  have  been  used  that  satisfy  some  or  all  of 
the  requirements  of  the  processor  in  various  parametric  receiver 
configurations.  The  common  element  in  each  of  these  techniques 
is  their  dependence  upon  the  resonance  properties  of  quartz 
crystals.  In  two  of  the  methods,  crystal  bandpass  or  crystal  band 
elimination  (notch)  filters  are  used  to  suppress  or  reject  the 
high  amplitude  pump  frequency  component  (carrier)  while  simul¬ 
taneously  passing  the  modulation  sidebands.  In  the  third  tech¬ 
nique,  the  modulation  sidebands  are  recovered  by  mixing  the 
parametric  receiver  signal  with  a  reference  signal  from  a  crystal 
stabilized,  phase  locked  loop  oscillator  in  a  high  performance 
diode  mixer.  Each  of  these  techniques  has  advantages  and  dis¬ 
advantages  that  will  be  discussed  in  the  following. 


3.2.1.  Bandpass  receiver.  It  is  well  known  that  the  infor¬ 
mation  in  a  modulated  carrier  is  redundantly  represented  in  the 
two  sidebands.  Single  sideband  (SSB)  communication  systems  have 
used  this  fact  to  advantage  for  many  years.  Thus  in  the  case  of 
the  parametric  receiver,  it  is  sufficient  to  recover  only  one 
sideband  signal.  One  method  of  accomplishing  this  is  to  use  a 
SSB  crystal  filter  to  reject  the  carrier  and  one  sideband  while 
passing  the  other  sideband.  Another  method  is  to  use  bandpass 
crystal  filters  that  pass  only  that  portion  of  one  sideband  that 
is  of  interest. 


The  latter  technique  was  used  in  early  experiments  to 
investigate  the  self-noise  sources  in  the  parametric  receiver  (17). 
A  block  diagram  of  the  receiver  used  in  those  experiments  is  shown 
in  Fig.  12.  The  skirts  of  the  cascaded  250  Hz  bandpass  crystal 
filters  suppressed  the  high  level  carrier  and  the  low  noise  ampli¬ 
fiers  provided  buffering  and  sufficient  amplification  of  the  band- 
passed  SSB  signal  to  drive  a  balanced  modulator.  The  demodulated 
signal  was  low  pass  filtered  to  yield  the  baseband  signal.  This 
receiver  had  a  5  dB  noise  figure  at  an  input  impedance  of  10  kft 
and  was  used  with  a  60  kHz  crystal  controlled  pump  signal  to 
receive  acoustic  signals  over  the  125  to  375  Hz  frequency  range. 

I 

The  bandpass  receiver  is  relatively  simple  and  easy  to 
implement  because  low  noise  integrated  circuit  preamplifiers  and 
wideband  balanced  modulators  are  readily  available.  A  minimum 
of  electronic  circuitry  is  required  to  translate  the  SSB  signal 
to  baseband.  The  bandpass  receiver  is  useful  for  some  experi¬ 
mentation  and  perhaps  for  limited  applications;  however,  it  has 
several  significant  disadvantages. 

One  disadvantage  of  the  bandpass  receiver  technique  is  that 
the  shape  factor  of  SSB  or  bandpass  crystal  filters  limits  the 
carrier  suppression  that  can  be  obtained  when  the  passband  ex¬ 
tends  to  within  a  few  hertz  of  the  carrier  frequency.  This  limits 
the  low  frequency  response  of  the  receiver  and,  hence,  the  para¬ 
metric  receiver.  Another  disadvantage  is  the  maximum  carrier-to- 
spectrum  sideband  noise  ratio  (C/SN)  that  can  be  achieved.  In 
the  bandpass  receiver  described  previously,  the  maximum  C/SN  ob¬ 
served  during  experiments  was  152  dB.  The  bandpass  receiver  was 
also  very  sensitive  to  vibration  induced  spurious  noise  at  high 
carrier  levels. 

3.2.2.  Band  elimination  receiver.  Band  elimination  crystal 
filters  were  first  used  for  carrier  suppression  in  parametric 
receiver  experiments  by  Barnard  et  al.  (22).  The  amplitudes  of 
high  level  pulsed  sinusoids  used  in  these  experiments  were  detected 
to  plot  beam  patterns  or  to  compare  with  predicted  interaction 
component  amplitudes.  The  sidebands  were  not  demodulated  to  yield 
a  broadband  receiver  and  signal  frequencies  below  1000  Hz  were 
not  investigated. 

The  band  elimination  filter  approach  was  combined  with  con¬ 
cepts  from  the  bandpass  receiver  in  developing  a  band  elimination 
receiver  capable  of  detecting  sideband  signals  over  the  35  to 
4000  Hz  frequency  range  at  carrier-to-sideband  ratios  approaching 
180  dB.  Figure  13  shows  a  block  diagram  of  that  receiver,  which 
performs  the  basic  functions  of  carrier  suppression,  sideband 
signal  amplification,  frequency  translation,  and  sideband  separa¬ 
tion.  The  band  elimination  receiver  is  functionally  similar  to 
the  bandpass  receiver  described  previously  but  the  change  to 
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Figure  13-  Block  diagram  of  the  band  elimination  receiver. 


band  elimination  filters  required  a  more  complex  demodulation 
process  to  accommodate  the  double  sideband  (DSB)  signal. 

Carrier  suppression  and  sideband  amplification  are  provided 
by  cascaded  crystal  filters  and  low  noise  amplifiers.  Each  filter 
provides  greater  than  40  dB  of  attenuation  at  frequencies  less  than 
10  Hz  from  the  65  kHz  carrier  frequency  while  signals  with  fre¬ 
quencies  greater  than  35  Hz  from  the  carrier  frequency  are  attenu¬ 
ated  less  than  3  dB.  Insertion  loss  of  each  crystal  filter  is  less 
than  2  dB  over  the  passband.  The  low  noise  amplifier  provides 
amplification  with  less  than  1  dB  of  additive  noise  to  the  signal 
from  the  crystal  filter  terminated  by  a  15  kft  resistor.  This 
combination  gives  a  receiver  noise  figure  of  5  dB  for  the  15  kft 
input  impedance  of  the  receiver. 

The  cascaded  crystal  filters  attenuate  the  carrier  by  greater 
than  140  dB  to  yield  a  double  sideband  suppressed  carrier  signal  at 
the  output  of  the  last  crystal  filter.  To  separate  the  sideband 
signals,  this  suppressed  carrier  signal  is  amplified  and  coupled 
into  a  pair  of  balanced  modulators  operating  in  phase  quadrature. 
The  signals  from  the  modulators  are  low  pass  filtered  to  remove 
the  higher  order  product  frequencies  and  any  residual  carrier. 

The  low  passed  signals  are  phase  shifted  by  a  broadband  quadrature 
phase  shift  of  90°  between  the  two  signal  channels.  The  signals 
from  the  phase  shift  networks  are  combined  with  appropriate  polari¬ 
ties  to  yield  the  upper  and  lower  sideband  signals  simultaneously. 

A  phase  locked  loop  (PLL)  oscillator  generates  the  carrier 
frequency  quadrature  reference  signals  for  the  balanced  modulators. 
The  PLL  oscillator  contains  a  digital  divider  network  and  a  volt¬ 
age  controlled  oscillator  operating  at  four  times  the  carrier  fre¬ 
quency.  The  final  stages  of  the  divider  network  are  driven  with 
complementary  waveforms  to  produce  the  carrier  frequency  quadrature 
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reference  signals.  The  carrier  frequency  reference  signal  for  the 
PLL  oscillator  is  normally  supplied  by  an  internal  crystal  stabi¬ 
lized  local  oscillator;  however,  a  switch  selectable  input  allows 
an  optional  external  reference  signal  to  be  used  if  needed. 

The  equivalent  input  noise  of  the  band  elimination  receiver 
was  measured  at  -151  dB  re  1  V//Hz  over  the  35  to  4000  Hz  fre¬ 
quency  range.  Maximum  power  dissipation  in  the  crystal  filter  is 
1  mW,  which  corresponds  to  a  maximum  carrier  level  of  29  dB  re 
1  V  from  a  15  kft  source.  Thus  the  receiver  should  be  able  to  de¬ 
tect  signals  180  dB  below  the  maximum  carrier  level.  During 
development  and  testing  of  a  340  m  parametric  receiver  at  Lake 
Travis  (Texas)  signals  as  small  as  -168  dB  referenced  to  the 
carrier  level  were  received  and  displayed  (23).  Rohde  et  al. 
describe  this  receiver  in  greater  detail  and  give  test  results 
from  the  development  (24) . 

The  band  elimination  receiver  discussed  in  this  section  is  a 
very  effective  processor  and  it  satisfies  the  requirements  of  the 
parametric  receiver  under  a  variety  of  conditions.  It  is  capable 
of  detecting  signals  with  smaller  modulation  indices  than  any 
other  currently  available  receiver.  The  simultaneous  availability 
of  both  sideband  signals  is  especially  useful  in  some  investiga¬ 
tions  (25).  The  most  serious  disadvantages  of  this  receiver  are 
the  high  input  impedance  of  the  crystal  filters  and  sensitivity  of 
these  filters  to  spurious  noise  induced  by  vibration  of  the  filters 
at  high  carrier  levels.  The  linear  characteristics  of  the  band 
elimination  receiver  can  be  a  disadvantage  in  applications  where 
the  transducers  are  subjected  to  large  amplitude  vibrations.  This 
effect  and  a  technique  for  reducing  the  detrimental  effects  of 
transducer  vibration  are  discussed  in  the  next  section. 

3.2.3.  Phase  locked  loop  receiver.  Because  the  water  is 
only  weakly  nonlinear,  in  virtually  all  cases  of  practical  in¬ 
terest  the  acoustic  phase  modulation  index,  <|>a(t),  is  very  small 
and  Eq.  (13)  can  be  simplified  to 

s ( t )  =  cos(u>ot+t|0  -  <f>a(t)  sin(wot+if0  .  (20) 

Reeves  et  al.  (26)  have  extended  the  analysis  to  include  trans¬ 
ducer  vibration.  If  the  pump  and  hydrophone  are  not  stationary 
relative  to  each  other,  the  phase  shift  tj)  is  not  a  constant  but 
is  given  by 

<Kt)  =  +  4>m(t)  >  (21) 

where  t pQ  is  the  phase  shift  due  to  propagation  delay  between  the 
mean  positions  of  the  pump  and  hydrophone,  and  <f>m(t)  is  the  phase 
variation  caused  by  displacements  of  the  transducers  from  their 
mean  positions.  In  contrast  to  |<f>a(t)|,  which  is  small  in  almost 


all  cases  of  interest,  |4>m(t)|  may  or  may  not  be  small  in  many 
applications,  depending  upon  the  vibration  amplitudes. 

If  |<t>m(t)|«l,  Eq.  (13)  can  be  simplified  as  previously 
to  yield 

s ( t )  =  cos (uot+0o)  -  [<i>a(t)+<j>m(t)  ]  sin (aJ0t+ii>0)  .  (22) 

In  this  case  the  modulating  signal  is  the  sum  of  the  acoustic 
modulation  due  to  nonlinearity  of  the  water  and  the  vibration 
modulation  due  to  transducer  motion.  Under  these  conditions  a 
linear  receiver  such  as  the  band  elimination  receiver  is  adequate. 

On  the  other  hand,  if  |<j>m(t)|  is  not  small  relative  to  a 
radian,  Eq.  (13)  can  only  be  linearized  in  (J>  (t),  which  gives 

Si 

s(t)  =  cos(0Jot+4>m(t)+iJ;o)  -  <f>a(t)  sin(a)0t+<j>m(t)+ij>0)  .  (23) 

Reeves  et  al.  (26)  have  shown  that  large  amplitude  transducer 
vibration,  which  produces  a  convolution  of  the  vibration  modulated 
carrier  spectrum  with  the  acoustic  spectrum  as  shown  by  Eq.  (23), 
causes  distortion  of  the  acoustic  signal  when  received  by  a  linear 
receiver  such  as  the  band  elimination  receiver.  A  phase  locked 
loop  technique  was  suggested  as  a  possible  solution  to  this 
problem. 

Investigation  of  the  phase  locked  loop  receiver  technique 
has  been  continued  by  Lamb  (27).  He  has  analyzed  the  carrier 
tracking  phase  locked  loop  and  designed  a  PLL  demodulator  for  use 
with  the  parametric  receiver.  A  block  diagram  of  this  PLL  re¬ 
ceiver  is  shown  in  Fig.  14.  The  PLL  receiver  operates  as  a 
closed  loop  servo  system  with  signal  phase  as  the  controlled 
variable.  The  signal  from  the  hydrophone  is  impedance  trans¬ 
formed  and  coupled  into  the  double  balanced  diode  mixer  which  has 
a  very  large  dynamic  range.  The  output  of  the  mixer  is  amplified 
and  low  pass  filtered  to  provide  the  control  voltage  for  the  low 
noise  voltage  ■"ontrolled  oscillator,  which  in  turn  provides  the 
carrier  frequency  reference  signal  for  the  mixer.  The  accuracy 
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Figure  14.  Block  diagram  of  the  phase  locked  loop  receiver. 
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with  which  the  voltage  controlled  oscillator  tracks  the  phase 
variations  of  the  input  signal  is  a  function  of  the  loop  gain  and 
the  variable  low  pass  filter  allows  the  loop  bandwidth  to  be 
changed  easily.  Tests  with  a  crystal  stabilized  low  noise  oscil¬ 
lator  showed  that  the  PLL  receiver  noise  floor  was  -160  dB 
referenced  to  the  carrier  level. 

The  performance  of  the  PLL  receiver  was  compared  to  that  of 
the  band  elimination  receiver  in  processing  the  output  of  a  para¬ 
metric  receiver  hydrophone  undergoing  large  amplitude,  low  fre¬ 
quency  vibration.  The  output  of  the  band  elimination  receiver  is 
shown  in  Fig.  15(a).  A  70  Hz  acoustic  signal  was  being  received 
while  the  parametric  receiver  hydrophone  was  mechanically  vibrated 
at  a  frequency  of  0.53  Hz  with  a  peak  displacement  of  4  mm.  The 
vibration  induced  modulation  sidebands  about  the  70  Hz  signal  are 
clearly  evident.  The  output  of  the  PLL  receiver  processing  the 
signal  from  the  parametric  receiver  under  the  same  conditions  is 
shown  in  Fig.  15(b).  The  improvement  in  the  intermodulation  side¬ 
bands  is  approximately  20  dB  for  these  data  which  were  obtained 
with  a  10  Hz  loop  filter.  Greater  reductions  of  the  intermodula¬ 
tion  sidebands  were  obtained  with  a  narrower  loop  filter.  Clearly 
under  these  conditions,  detection  of  narrowband  signals  would  be 
more  easily  accomplished  with  the  PLL  receiver. 

The  major  disadvantage  of  the  PLL  receiver  at  this  time  is 
the  maximum  carrier-to-sideband  ratio  of  160  dB.  An  improvement 
of  10  to  20  dB  in  this  parameter  is  needed  for  some  applications. 
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Figure.  15.  Comparison  of  performance  of  band  elimination 
and  phase  locked  loop  receivers  with  large  amplitude  trans¬ 
ducer  vibration.  A  70  Hz  acoustic  signal  was  being  received 
while  the  hydrophone  was  vibrating  at  a  frequency  of  0.53  Hz 
and  a  peak  displacement  of  4  mm. 
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CAPT  Keathley,  Code  OP-954 
CAPT  W.  Christensen,  Code  OP-951F 
CAPT  R.  B.  Gilchrist,  OP-95T 
CAPT  J.  Van  Metre,  Code  OP-224 
LCDR  J.  Carter,  Code  OP-961C4 
CDR  Harlett ,  Code  OP-952D1 
CDR  Kelley 
Ms.  J.  C.  Bertrand 

Office  of  the  Chief  of  Naval  Operations 
Long  Range  Planning  Group 
2000  North  Beauregard  St. 

Alexandria,  VA  22311 

Attn:  CAPT  J.  R.  Seesholtz,  Code  00X1 
Director 

Naval  Research  Laboratory 
Department  of  the  Navy 
Washington,  DC  20375 
Attn:  J.  Munson,  Code  5100 

F.  Ingenito,  Code  5120 
W.  Newbauer,  Code  5108 

L.  Dragonette 

M.  Potosky,  Code  5109 
J.  Jarznski,  Code  5131 
R.  D.  Corsaro,  Code  5131 
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Commander 

Naval  Ship  Weapon  Systems  Engineering  Station 
Department  of  the  Navy 
Port  Hueneme,  CA  93043 
Attn:  Dr.  G.  Hetland,  Jr. 

Naval  Research  Laboratory 
Underwater  Sound  Reference  Division 
P.  0.  Box  8337 
Orlando,  FL  32856 
Attn:  J.  E.  Blue 

A.  L.  Van  Buren 
P.  H.  Rogers 


Commander 

Naval  Coastal  Systems  Center 
Department  of  the  Navy 
Panama  City,  FL  32407 
Attn:  H.  Warner,  Code  791 
D.  Skinner,  Code  791 
M.  Lacey,  Code  720 
J.  Hammond,  Code  722 
D.  Folds 

Officer-in-Charge 

New  London  Laboratory 

Naval  Underwater  Systems  Center 

Department  of  the  Navy 

New  London,  CT  06320 

Attn:  M.  Moffett 

R.  H.  Mellen 
W.  Konrad 
R.  Schumaker 
L.  King 

Commanding  Officer 
Naval  Surface  Weapons  Center 
White  Oak  Laboratory 
Department  of  the  Navy 
Silver  Spring,  fID  20910 
Attn:  G.  Gaunaurd,  Code  R31 

G.  Kalaf,  Code  U5 

Commander 

Naval  Ocean  Systems  Center 
Department  of  the  Navy 
San  Diego,  CA  92152 
Attn:  H.  Bucker,  Code  5311 

F.  Gordon,  Code  633 

J.  Reeves,  Code  635 

W.  Angelof f ,  Code  6233 
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80 


Commanding  Officer 

Naval  Ocean  Research  and  Development  Activity 
NSTL  Station,  MS  39529 
Attn:  R.  Goodman,  Code  110 
S.  Marshal,  Code  340 
J.  Posey,  Code  340 
Mr.  Soloman,  Code  500 
A.  L.  Anderson,  Code  320 


Commander 

Naval  Intelligence  Command 
2461  Eisenhower  Avenue 
Alexandria,  VA  22331 

Commander 

Naval  Explosive  Ordnance  Disposal  Facility 
Indian  Head,  MD  20640 

Attn:  Mr.  L.  H.  Dickinson,  Technical  Director 

Commanding  Officer 

USCG  Research  and  Development  Center 

Avery  Point 

Groton,  CT  06340 

Attn:  CAPT  M.  Y.  Suzich 

Officer-in-Charge 

David  W.  Taylor  Naval  Ship  Research  and  Development  Center 
Department  of  the  Navy 
Bethesda,  MD  20084 
Attn:  R.  Eippeon 

National  Academy  of  Sciences 
Naval  Studies  Board 
2101  Constitution  Avenue,  NW 
Washington,  DC  20418 
Attn:  L.  Hunt 

Commanding  Officer 
Naval  Postgraduate  School 
Monterey,  CA  93940 
Attn:  Prof.  G.  Sacklman 

Prof.  A.  B.  Coppens 
Prof.  J.  V.  Sanders 
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North  Atlantic  Treaty  Organization 
SACLANT  ASW  Research  Centre 
APO  09019 
New  York,  NY 
Attn:  w.  Kuperman 
F.  Whicker 
F.  Jensen 

Applied  Research  Laboratory 
The  Pennsylvania  State  University 
Post  Office  Box  30 
State  College,  PA  16801 
Attn:  F.  H.  Fenlon 

The  Catholic  University  of  America 
Washington,  DC  20064 
Attn:  Dr.  Uberall 

Westinghouse  Electric  Corporation 
P.  O.  Box  1488 
Annapolis,  MD  21404 
Attn:  Mr.  A,  Nelkin 

Dr.  P,  J.  Welton 

Raytheon  Company 
P.  0.  Box  360 
Portsmouth,  RI  02871 
Attn:  Mr.  J.  F.  Bartram 

Mr.  E.  P.  Aurand 
19  Hanapepe  Place 
Honolulu,  HI  96825 

Department  of  Engineering  and  Applied  Science 

Yale  University 

New  Haven,  CT  06520 

Attn:  Dr.  P.  M.  Schultheiss 

RAMCOR,  Inc. 

800  Follin  Lane 
Vienna,  VA  22180 
Attn:  Mr,  V.  J.  Lujetic 

Systems  Planning  Corporation 
1500  Wilson  Blvd.,  Suite  1500 
Arlington,  VA  22209 
Attn:  Mr.  Jack  Fagan 
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102-113 


114 


Bolt,  Beranek  &  Newman,  Inc. 

50  Moulton  Street 
Cambridge,  MA  02138 
Attn:  Dr.  J.  E.  Barger 
Dr.  F.  J.  Jackson 

Scripps  Institution  of  Oceanography 
Marine  Physical  Laboratory 
University  of  California 
San  Diego,  CA  92152 
Attn:  Dr.  William  S.  Hodgkiss 

Tracor,  Inc. 

6500  Tracor  Lane 
Austin,  TX  78721 
Attn:  Mr.  J.  D.  Williams 
Mr .  J .  Dow 

Radian  Corporation 
8500  Shoal  Creek  Blvd. 

P.  0.  Box  9948 
Austin,  TX  78757 
Attn:  Mr.  Jerry  L.  Bardin 
Dr.  C.  R.  Reeves 

General  Physics  Corporation 
10630  Little  Patuxent  Parkway 
Columbia,  MD  21044 
Attn:  Dr.  Frank  Andrews 

Trans  World  Systems,  Inc. 

1311A  Dolly  Madison  Blvd. 

McLean,  VA  22101 
Attn:  Mr.  Sam  Francis 

Commanding  Officer  and  Director 
Defense  Technical  Information  Center 
Cameron  Station,  Building  5 
5010  Duke  Street 
Alexandria,  VA  22314 

Battelle  Memorial  Institute 
505  King  Avenue 
Columbus,  OH  43201 
Attn:  TACTEC 
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Office  of  Naval  Research 
Resident  Representative 
Room  582,  Federal  Building 
Austin,  TX  78701 

Physical  Acoustics  Group,  ARL:UT 

Garland  R.  Barnard,  ARL:UT 

C.  Robert  Culbertson,  ARL:UT 

Tommy  G.  Goldsberry,  ARL:UT 

Loyd  Hampton,  ARL:UT 

John  M.  Huckabay,  ARL:UT 

Robert  A.  Lamb,  ARL:UT 

T.  G.  Muir,  ARL : UT 

David  F.  Rohde,  ARL:UT 

Reuben  H.  Wallace,  ARL:UT 

Library,  ARL:UT 
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